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ABSTRACT 

We have carried out an extensive radio study with the Very Large Array on the 
Seyfert 1.5 galaxy Mrk 6 and imaged a spectacular radio structure in the source. 
The radio emission occurs on three different spatial scales, from ~7.5 kpc bubbles 
to ~1.5 kpc bubbles lying nearly orthogonal to them and a ~1 kpc radio jet lying 
orthogonal to the kpc-scale bubble. To explain the complex morphology, we first 
consider a scenario in which the radio structures are the result of superwinds 
ejected by a nuclear starburst. However, recent Spitzer observations of Mrk 6 
provide an upper limit to the star formation rate (SFR) of ~ 5.5 M yr _1 , an 
estimate much lower than the SFR of ~ 33 M & yr _1 derived assuming that the 
bubbles are a result of starburst winds energized by supernovae explosions. Thus, 
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a starburst alone cannot meet the energy requirements for the creation of the 
bubbles in Mrk 6. We show that a single plasmon model is energetically infeasible, 
and we argue that a jet-driven bubble model while energetically feasible does not 
produce the complex radio morphologies. Finally, we consider a model in which 
the complex radio structure is a result of an episodically-powered precessing jet 
that changes its orientation. This model is the most attractive as it can naturally 
explain the complex radio morphology, and is consistent with the energetics, the 
spectral index and the polarization structure. Radio emission in this scenario is 
a short-lived phenomenon in the lifetime of a Seyfert galaxy which results due to 
an accretion event. 

Subject headings: galaxies: Seyfert - - galaxies: individual (Mrk 6) - - radio 
continuum: galaxies 

1. INTRODUCTION 

Radio observations show that most Seyfert galaxies have sub-parsec-scale radio emission 
(e.g., de Bruyn & Wilson 1976; Ulvestad & Wilson 1984a,b; Roy et al. 1994; Thean et al. 
2000; Lai et al. 2004). Some of these galaxies with small-scale radio emission show elongated 
structures, similar to the radio jets seen in powerful radio galaxies. Radio emission on 
kiloparsec scales has also been detected in some Seyfert galaxies. The large-scale radio 
structures however are not typically found to be aligned with the small scale jets (Baum 
et al. 1993; Colbert et al. 1996; Gallimore et al. 2006). Further, Seyfert radio jets were found 
to be randomly aligned with respect to the host galaxy major axes (Kinney et al. 2000; 
Schmitt et al. 2001; Schmitt & Kinney 2002). Abrupt changes of jet axes were also found 
in some Seyfert galaxies, e.g., NGC 4151 (Ulvestad et al. 1998) and NGC 1068 (Gallimore 
et al. 1996). It appears that the radio structures in Seyfert galaxis are more complex than 
in the radio galaxies. Extensive discussions of large scale radio structures in Seyfert galaxies 
can be found in Baum et al. (1993); Colbert et al. (1996) and Gallimore et al. (2006). 

Mrk 6 is an early-type (SOa) Seyfert 1.5 galaxy. It is also one of the Seyfert galaxies that 
have radio emission on both small and large scales. Using the Westerbork Synthesis Radio 
telescope (WSRT), Baum et al. (1993) detected a pair of large weak "radio lobes" extending 
northeast-southwest and a smaller scale radio structure extending east- west. High resolution 
MERLIN observations of the nuclear region by Kukula et al. (1996) revealed a well-defined 
jet extending in the north-south direction. This jet is nearly aligned with the ionization cone 
reported earlier by Meaburn et al. (1989), and the optical polarization position angle (P.A.) 
derived through spectropolarimetric observations (Smith et al. 2004). Kukula et al. (1996) 
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also found some evidence of a pair of lobes on the sides of the jet which corresponded to the 
east-west radio structure observed by Baum et al. (1993) and Nagar et al. (1999). However, 
due to the relatively poor resolution of WSRT (~ 3". 5) and the limited u-v coverage of 
MERLIN, few conclusions could be made on the relation between the three scales of radio 
emission in this source. In order to have a better understanding of the radio structures 
in Mrk 6, we undertook a comprehensive radio study of Mrk 6 with the Very Large Array 
(VLA). 

Throughout this paper we assume Hq=71 km s _1 Mpc _1 , Q m = 0.27 and VL vac = 0.73. 
Therefore, at the redshift of z=0. 01881 (recession velocity ~5640 km sec" 1 ), or at the dis- 
tance of 80.6 Mpc for Mrk 6, 1 arcsec is equivalent to 377 parsecs. 



2. OBSERVATIONS AND DATA ANALYSIS 

Mrk 6 was observed with the VLA at A 6 cm (4860 MHz) and A 20 cm (1430 MHz) 
during 1995-1996. Almost all VLA configurations were used in the observations in order 
to image Mrk 6 on different spatial scales. 3C 286 was used as the primary flux density 
and polarization calibrator, while 0614+607 was used as the phase calibrator for the en- 
tire experiment. The data were processed with NRAO's software package AIPS using the 
standard imaging and self-calibration procedures. Table 1 lists the observing frequency, the 
corresponding bandwidth, the VLA configuration, the observation date, the FWHMs of the 
synthesized beams, the total flux densities of the source at various resolutions, and the r.m.s. 
noise in the final maps. 
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Table 1: Observation Log 



Frequency 


BW 


Config. 


Observation 


Resolution 


Total Flux Density 


r.m.s. 


(GHz) 


(MHz) 




Date 


(arcsec 2 ) 


(mJy) 


(^Jy/bm) 


4.86 


50 


A 


07/21/1995 


0.41x0.31 


85.8 


16 


4.86 


50 


B 


11/02/1995 


1.43x1.13 


96.4 


16 


4.86 


50 


C 


02/17/1996 


4.27x4.04 


96.5 


25 


4.86 


50 


D 


05/03/1995 


16.2x11.4 


100.5 


29 


1.43 


50 


A 


07/03/1995 


1.44x1.12 


274.4 


27 


1.43 


50 


C 


02/17/1996 


13.9x13.0 


270.8 


147 


1.43 


50 


D 


06/01/1995 


59.1x38.2 


275.0 


210 
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3. RESULTS 

The total intensity radio maps are presented in Fig. 1 while Fig. 2 shows the relative 
orientation of the radio structures with respect to the extended emission- line region (ENLR) 
and host galaxy of Mrk 6. Figures 3 and 4 show the radio polarization and the spectral 
index images, respectively. We describe below the primary findings of our study. 

3.1. Radio Morphology: Emission on Three Spatial Scales 

Figure 1 displays the mosaic of radio emission at 6 cm. The radio map of the inner jets is 
from the MERLIN observations of Kukula et al. (1996). It is immediately apparent that Mrk 
6 has radio emission on three scales. First, there is a radio jet of extent ~1 kpc (3") extending 
in the north-south direction, at the active galactic nucleus (AGN). Second, there are a pair 
of inner bubbles of size ~1.5 x 1.5 kpc (~ 4" x 4") each, extending east- west and nearly 
perpendicular to the central jets. And third, there are a pair of outer bubbles of dimensions 
~4.5 x 7.5 kpc (~ 12" x 20") each, extending northeast-southwest. In Fig. 2 we observe the 
relative orientation of these radio structures with respect to the extended [OIII] emission-line 
region as observed by Kukula et al. (1996) and the host galaxy optical continuum from the 
high-resolution HST image (obtained from the HST archive). We note that the inner kpc 
jet is roughly aligned with the emission-line gas. Further, the outer bubbles are aligned with 
the host galaxy minor axis. The SDSS image from the NASA/IPAC Extragalactic Database 
(NED) indicates that the extent of the galaxy is > 15 x 20 kpc; the radio structures therefore 
lie within the confines of the host galaxy. 
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Fig. 1. — Grey scale images of the 6 cm radio structures observed in Mrk 6. The bubble- 
like structures are clearly seen in the VLA images. The MERLIN map of the inner jet is 
taken from Kukula et al. (1996). The centre of the optical host galaxy lies close to the (0,0) 
position in the MERLIN image, with a 3a uncertainity of 0.27" (Clements 1981). There is 
no radio 'core' coincident with this position. 
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We have listed some geometric and physical parameters pertaining to these individual 
radio structures, including their sizes, position angles, flux density and spectral indices in 
Table 2. For the outer pair of bubbles, the radio flux density of the northeast bubble is 
about two and a half times that of the southwest bubble. For the inner pair, the radio flux 
density of the east bubble is about two times that of the west bubble. We note that (1) 
the two pairs of bubbles and the jets show self-similar characteristics in their morphologies 
and sizes, i.e., the width of the larger structure corresponds to the overall size of the next 
smaller structure, (2) the bubbles show clear edge-brightening and (3) the bubbles contain 
filamentary structures. 



- 8- 



1.0 1.5 2.0 2.5 




06 45 48 47 46 45 44 43 42 41 40 
RIGHT ASCENSION 

5 10 15 




06 52 16 15 14 13 12 11 10 09 

RIGHT ASCENSION 



Fig. 2. — The 20 cm radio contours superimposed on the grey-scale images of (top) extended 
[OIII] emission-line region (Kukulaet al. 1996), and (bottom) HST (WFPC2/F606W) optical 
continuum from the host galaxy (obtained from the HST archive). The outer bubble-like 
structures are projected roughly orthogonal to the major axis of the bulge-dominated SOa 
host galaxy. However, the structures lie within the extent of the host galaxy. 
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Table 2: Radio Properties of Mrk 6 



Structure 


Size 

(arcsec 2 ) 


P.A. 


s 6 

(mJy) 


L 6 

(W Hz" 1 ) 


af 


Jet 


1x3 


179° 


75.0 


5.8E22 


-0.91±0.21 


North 










-0.90±0.14 


South 










-0.84±0.10 


Inner Bubble 




76° 


11.4 


8.9E21 




East 


4x4 




7.9 


6.1E21 


-0.61±0.37 


West 


3x3 




3.5 


2.7E21 


-0.55±0.30 


Outer Bubble 




33° 


9.4 


7.3E21 




North 


12x20 




6.2 


4.8E21 


-0.71±0.21 


South 


12x18 




3.2 


2.5E21 


-0.55±0.26 



Note. — Col. 2: Extents of structures using the 6 cm maps; Col. 4: 6 cm flux density in mJy; Col. 5: 6 cm 
luminosity in W Hz -1 ; Col. 6: Spectral index estimated using the 6 cm and 20 cm images. 



3.2. Radio Polarization 



Polarization was detected at 6 cm in the region corresponding to the western edge of 
the inner bubble and on the edges of the outer pair of bubbles. Figure 3 displays the 6 cm 
VLA C-configuration image with polarization electric vectors superimposed. The length of 
the vectors is proportional to the fractional polarization. At this configuration, polarization 
was detected at the 10-20 a levels (or 0.2-0.4 mJy beam -1 ). This corresponds to a high 
fractional polarization of ~ 50% at the edges and less than 1% polarization near the center. 

We were unable to detect polarization on the edge of the outer bubbles with the B- 
configuration, and anywhere with A-configuration, due to the smaller beam size and lower 
resultant surface brightness sensitivity. A suitable taper on the B-configuration -ut>-coverage 
recovers the polarization features observed in the C-configuration image (Fig. 3), thus con- 
firming the detection of polarization. Assuming that the radio emitting medium is optically 
thin, and the Faraday rotation is small at 6 cm (Galactic rotation measure is < 30 rad m -2 
in the direction of Mrk 6, Simard-Normandin & Kronberg 1980), we infer the magnetic field 
orientation to be mostly aligned with the edge of the bubbles. 



20 15 10 5 -5 -10 -15 -20 

ARC SEC 



Fig. 3. — Contour image of Mrk 6 at 6 cm with polarization electric vectors and percentage 
of polarization indicated. The image was made using the VLA C-configuration data. The 
contour levels are 7.45 l(r 5 x (-1,1,2,4,6,8,12,16,24,32,64,128, 256,512) Jy/beam, with the 
beam size 4". 2 x 4".0. The lowest contour level corresponds to 3 times the off source r.m.s 
noise. 1" in length of polarization vector corresponds to 10% of polarization. 
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3.3. Spectral Indices 

Figure 4 (top) shows the 6 cm B-configuration contour image overlaid on the grey scale 
spectral index image, which was made from the 6 cm B-configuration and the 20 cm A- 
configuration images, after restoring both images with the same synthesized beam. The 
spectral indices vary only slightly across the image. The average spectral index over the 
central region corresponding to the inner bubbles is — 0.78±0.26, that over the northern 
bubble is — 0.71±0.21 and that over the southern bubble is — 0.55±0.26. These spectral 
indices are consistent given the large errors. 

Figure 4 (bottom) displays the spectral index map between 6 cm and 18 cm, in the 
neighbourhood of the jets. This image was made by combining our 6 cm A-configuration 
map with the 18 cm MERLIN map from Kukula et al. (1996), at the matched resolution of 
0".4 x 0".4. It shows that the spectral index of the northern hotspot in the jet is — 0.90±0.14, 
that of the southern hotspot is — 0.84±0.10 and that of the saddle between the two hotspots 
is — 0.96±0.08. Our results are in general agreement with those of Kukula et al. (1996). The 
average spectral index over the eastern bubble is — 0.61±0.37 and over the western bubble 
is — 0.55±0.30, slightly flatter than those in the jets, but similar to those in the outer pair 
of bubbles. 



3.4. Energetics 

We have estimated the magnetic field strengths in regions corresponding to the bubbles 
assuming equipartition of energy between relativistic particles and the magnetic field (Bur- 
bidge 1959). Using Eqns. 1-5 of O'Dea & Owen (1987), we have also obtained the minimum 
pressure and particle energy (electrons and protons) at minimum pressure. These estimates 
are listed in Table 3. The averaged total radio luminosity (L rad ) was estimated assuming 
that the radio spectrum extends from 10 MHz to 10 GHz with a spectral index of a=-0.6. 
Further, it was assumed that the relativistic protons and electrons have equal energies. 

In Table 3 we list the minimum pressure values for the two scenarios where the bubbles 
are completely filled (volume filling factor, 0=1) or are mostly empty (0=1O~ 3 ). The esti- 
mated minimum pressure in the filled outer bubbles turns out to be ~ 10~ 12 dynes cm~ 2 and 
in the inner bubbles P m j„ is ~ 10~ n dynes cm -2 . The pressure is about fifty times larger 
for the bubbles when the volume filling factor is 0=1O~ 3 . 
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Fig. 4. — Spectral index map (see text) in greyscale overlaid on the (top) 20 cm A- 
configuration and (bottom) 6 cm A-configuration contour maps. The contour levels 
and beamsizes are (top) 8.82 10" 5 x (-1,1,2,4,6,8,12,16,24,32,64,128,256,512) Jy/beam and 
4".5 x 4".5 and (bottom) 7.0 10~ 5 x (-1,1,2,4,6,8,12,16,24,32,64,128,256,512) Jy/beam and 
0".5 x 0".5. 
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4. DISCUSSION 

Bubbles and rings have been observed in a few radio galaxies and quasars, for example, 
in Hercules A (Dreher & Feigelson 1984; Baum et al. 1996), 3C 310 (van Breugel & Fomalont 
1984), and Centaurus A (Quillen et al. 2006). Bubble-like structures have also been observed 
in some Seyfert and starburst galaxies, e.g., NGC 2992, NGC 3079, NGC 5548, NGC 4051, 
and NGC 6764 (Wehrle & Morris 1988; Duric & Seaquist 1988; Baum et al. 1993; Nagar et al. 
1999), and planetary nebulae such as rj Carina and Hubble 5 (Nota et al. 1995; Dwarkadas & 
Balick 1998). Several mechanisms have been suggested for the formation of these bubble-like 
structures. We describe them in Sect. 4.2 and subsequently examine three potentially viable 
models for the formation of the radio structures in Mrk 6. 



4.1. Constraints On A Model For Mrk 6 

We note that the bubble-like structures observed in the above mentioned radio galaxies 
are essentially spherical and they are either part of or directly connected with the radio jet. 
This is different from the bubbles observed in Mrk 6 which are elongated, larger than the jet 
and not aligned with it. The elongated morphology of the bubbles make them more similar 
to the bubbles observed in planetary nebulae. 

In order to estimate the originating point of the bubbles (where the two edges meet, see 
Fig. 5), under the simplifying assumption that the two edges of a bubble on either side lie 
in the same plane, we used a simple IDL program which, given the pixel positions of certain 
bright spots on the edges of the bubbles, fits a third degree polynomial to them, and derives 
the point where the difference in the position of the fitted functions is a minimum. Thus, 
it turns out that the outer pair of bubbles originate near a=06:52: 12.359, <5=+74:25:37.26 
(J2000) while the inner pair of bubbles originate near a=06:52: 12.340, <5=+74:25:37.081 
(J2000). Both originate from near the same place (offset by 0".3 or ~110 pc), which inter- 
estingly, lies around the hot spots 3 and 4 in Kukula et al. (1996) (to the south of the (0,0) 
position in the MERLIN image of Fig. 1). It should be pointed out that the bright spots we 
pick on the edges of the bubbles are not without ambiguities, which in turn give uncertain 
positions for the expected origins. The offset in the originating point of the bubbles may 
therefore not be of much significance. As Kukula et al. (1996) have noted, no radio 'core' 
has been found coincident with the centre of the optical host galaxy, which lies close to the 
(0,0) position (with a 3a uncertainity in position of 0.27", Clements (1981)) in the MERLIN 
image (Fig. 1), and it is not clear if one of the knots ~ 0.5" to the south, is the radio core. 
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Fig. 5. — Extrapolations of bubbles down to their origins (see text). The axes coordinates 
are in pixels, the pixel size being 0.1" in the top and 0.25" in the bottom panel, respectively. 
This figure indicates that both pairs of bubbles originate from a region within 170 pc. 
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Table 3: Physical Properties and Energetics 



Structure, Filling 


Length 


Width 








p . 

1 mm 


r actor (0) 


///\ 

( ) 


( ) 


/ — 1 \ 

(erg s J 


(erg) 


(jiG) 


(dynes 


cm 2 ) 


Outer Bubble, 0=1 
















North 


20.4 


12.1 


7.4E38 


9.7E54 


3.6 


1.2E- 


•12 


South 


17.8 


12.4 


3.8E38 


6.4E54 


3.0 


8.7E- 


13 


Inner Bubble, 0=1 
















East 


3.8 


3.7 


9.4e38 


2.0E54 


12.3 


1.4E- 


11 


West 


3.0 


3.2 


4.2e38 


9.9E53 


11.3 


1.2E- 


11 


Outer Bubble, 0=10 6 
















North 






7.4E38 


5.1E53 


25.9 


6.2E- 


11 


South 






3.8E38 


3.3E53 


22.0 


4.5E- 


11 


T T — X 111 / -1 f\ X 

Inner Bubble, 0=10 
















East 






9.4E38 


1.0E53 


88.4 


7.2E- 


•10 


West 






4.ZilvOO 


O.lrL/OZ 


el. 4 


6.1E- 


•10 


North-South Jet 


3.1 


1.0 


9.0E39 


2.1E54 


52.3 


2.5E- 


•10 


Outer Jet (North) 


6.8 


3.0 


1.3E38 


6.9E53 


6.7 


4.1E- 


12 




7.5 


3.1 


1.5E38 


8.1E53 


6.7 


4.1E- 


12 


Outer Jet (South) 


7.2 


3.1 


8.7E37 


5.7E53 


5.7 


3.0E- 


12 




5.4 


2.0 


3.6E37 


2.1E53 


6.2 


3.6E- 


12 


Inner Jet (East) 


1.8 


0.8 


3.5E38 


2.2E53 


27.4 


6.9E- 


11 




1.5 


0.7 


7.3E37 


7.4E52 


19.9 


3.7E- 


•11 


Inner Jet (West) 


1.4 


0.6 


3.2E37 


3.9E52 


17.6 


2.9E- 


11 




1.5 


0.7 


1.8E38 


1.2E53 


25.8 


6.2E- 


11 



Note. — Col. 2 and Col. 3: The length and width in arcsecs of the volume element considered for the 
calculations, using; the 6 cm image. Lrad — averaged total radio luminosity (see Sect. 3.4), E m j n and B m j n 
correspond to the particle energy and magnetic field strength at the minimum pressure P m in, obtained 
assuming the 'equipartition' condition. The upper part of the table has the parameters listed for the bubbles 
which are either completely filled (0=1) or mostly empty (0=1O~ 3 ), while the lower part of the table lists the 
parameters for the radio jet and volume elements taken along the edges of the bubbles, assuming 0=1. In 
the precessing jet scenario, these volume elements correspond to portions of the precessing radio synchrotron 
jet. Outer Jet (North) corresponds to a volume element taken along the edge of the northern outer bubble 
on either side, and so on. 
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4.2. Mechanisms for Producing Bubble-like Structures 

Extended radio structures have been observed in both Seyfert and starburst galaxies, 
and some of the structures are bubble-like (Colbert et al. 1996; Heckman et al. 1993). In 
starburst galaxies, the bubble-like radio structures have been suggested to form through the 
interaction between the starburst superwind and the ISM (Duric & Seaquist 1988; Heckman 
et al. 1993). In the case of Seyfert galaxies, Colbert et al. (1996) found that the kpc-scale 
radio outflows are more consistent with diverted AGN-driven jets. But they didn't rule out 
the possibility that the outflows are starburst superwinds. In fact, indications of both AGN 
and starburst activity have been found in Seyfert galaxies like NGC 7469 (Genzel et al. 
1995), NGC 3079 (Veilleux et al. 1994) and the Circinus galaxy (Maiolino et al. 1998). 

Other scenarios have been proposed to explain bubble formation in galaxies. In order to 
explain the spherical shell-like structures in the radio galaxies Her A and 3C 310, Morrison 
& Sadun (1996) introduced a model wherein acoustic waves or weak shocks are excited along 
the jet axis of the AGN, in a pre-existing thermal galactic wind. These acoustic shells expand 
uniformly at the speed of sound and drift along with the wind. At a much later second stage, 
a new and much faster flow of relativistic plasma is energized by the AGN, which fills these 
shells with radio-emitting electrons, thus creating radio-bright bubbles. Alternately, Pedlar 
et al. (1985) invoked a model of expanding radio-emitting plasmons which shock, compress 
and accelerate the forbidden-line gas in Seyfert galaxies. Both these models require the 
bubbles to be continuously injected with relativitic particles from the AGN. 

Bubble-like structures observed in some planetary nebulae like Hb 5, (Terzian & Hajian 
2000) 1 or r] Carina (Nota et al. 1995), have been suggested to arise due to the presence of 
orbiting binary companions (Soker 2004). A binary black hole system could similarly be 
invoked to explain the bubble-like structures in Mrk 6. In this picture, the primary black 
hole is spinning and is responsible for the formation of the radio jet. The secondary black 
hole on the other hand, is not spinning but orbits the primary black hole in a plane which 
is perpendicular to the inner bubbles. The quasi-sherical disk wind from the primary black 
hole forms an accretion disk around the secondary black hole. This in turn launches fast 
winds into the pre-existing "density contrast" between the equatorial and polar directions 
(see Balick 1987) created by the orbiting secondary black hole. 

The complexity involved in this binary blackhole model with the two backholes hav- 
ing perpendicular orbital planes, however, makes it unattractive. For example, the orbital 
plane of the secondary black hole will be aligned with the inner jet axis and must exert a 
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disturbing influence on its propagation (Merritt & Milosavljevic 2005), an effect which is 
not observed in the high resolution MERLIN image. Further, the existence of a secondary 
black hole orbiting in a plane perpendicular to the accretion disk of the primary, must have 
discernable effects (Sillanpaa et al. 1988; Lehto & Valtonen 1996) on the dynamics of the 100 
pc-scale disk (Gallimore et al. 1998) which we assume to be relatively stable according to the 
ENLR (Meaburn et al. 1989, see Fig. 2), and the optical polarization position angle derived 
through spectropolarimetric observations (Smith et al. 2004). It is also uncertain whether 
the planetary nebula model will scale to an AGN. In the following sections, we discuss in 
some detail the pros and cons of three potentially viable scenarios that could give rise to the 
bubble-like structures in Mrk 6, viz., (i) a starburst-driven bubble scenario, (ii) a jet-driven 
bubble scenario, and (iii) a precessing jet scenario. 



4-2.1. Starburst-driven Bubble Scenario 

The morphology of the pairs of the bubbles, especially their edge-brightening suggests 
that the radio emission could come from the shells of superbubbles. The sharp and well- 
defined edges of the bubbles and the existence of filaments in the bubbles would be consistent 
with strong interaction between the bubbles and the interstellar medium (ISM) where the 
radio-emitting particles (relativistic electrons) are accelerated and/or the magnetic field gets 
amplified (Pacholczyk & Scott 1976). For example, if the expansion of the bubbles are 
caused by diffusion with little interaction (shocks), we would expect the bubbles to be more 
irregular and diffusive, as seen in most radio lobes. The spectral index of the bubbles seems 
rather flat, also consistent with interaction and particle re-acceleration. 

Assuming that the outer bubbles are formed through starburst winds energized by 
supernovae explosions, we estimate the energy budget. If we assume the bubbles are filled 
with the ejecta from the supernovae, the energy required for the bubbles to expand to the 
present sizes (V~ 7.5 x 4.5 x 4.5 kpc 3 ) adiabatically within the ISM of pressure P~ 10~ 10 
dynes cm -2 is 4PV=1.7xl0 57 ergs. This corresponds to 1.7 x 10 6 supernova explosions with 
individual explosive energy of 10 51 ergs. For wind speeds of 1000 km s^ 1 , the time required 
for the bubbles to expand to the present size is ~ 7 x 10 6 yrs. Thus the required supernova 
explosion rate is ~ 0.24 yr _1 , which is comparable to the predicted rate for the prototypical 
starburst galaxy M 82 (Seaquist & Odegard 1991; Rieke et al. 1993). The required supernova 
rate would be reduced if the bubbles are not fully filled with the superwinds. We must 
mention that we have not included the loss through radio emission since it is negligible 
compared with the dynamical energy we have estimated above. 

The derived supernova rate can yield an estimate of the required star-formation rate 
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(SFR). By assuming a Salpeter initial mass function over the mass range 0.1 < M/M ^ 100, 
using the above derived supernova rate, we obtain (by using Eqns. 28, Condon et al. (2002) 
and 20, Condon (1992)), a large star-formation rate of ~ 33 M yr -1 for Mrk 6. The SFR 
estimated from the 1.4 GHz luminosity of the bubbles (~ 6 x 10 22 W Hz -1 ) for a = —0.6 
and Eqns. 21 of Condon (1992) and 28 of Condon et al. (2002), turns out to be much larger 
(SFR -76 M yr- 1 ). 

An upper limit to the SFR can be obtained from recent Spitzer IRAC observations 
of Mrk 6 (Buchanan et al. 2006). Using the 8/xm flux density of 0.1 Jy resulting from 
the central ~ 7" region (corresponding to a spatial scale of 2.7 kpc), and the relation, 
SFR = 2.0xl0~ 43 v L u (8 \xm) erg sec -1 , (Calzetti et al. 2006, in prep.), which makes use 
of the nebular lines calibration of Kennicutt (1998), we obtain a SFR (M ^ O.1M ) ~ 5.5 
M yr -1 . This value is much lower than the SFR derived assuming that the bubbles are a 
result of starburst winds energized by supernovae explosions. 

The derived minimum pressure estimate for the bubbles (~ 10~ n — 10~ 12 dynes cm -2 ) 
is one to two orders of magnitude lower than the typical ISM pressure of ~ 10~ 10 dynes cm -2 
found in an elliptical or an SO galaxy (Mathews & Brighenti 2003). This could imply that 
either (i) the pressure has a large contribution from thermal particles, (ii) the equipartition 
condition does not exist in the bubbles, or (iii) the equipartition condition exists and the 
derived pressure estimates are the actual ISM pressures in Mrk 6. Further, through optical 
spectroscopic observations, Heckman et al. (1990) have found typical kpc-scale pressures of 
2 — 4 x 10 -9 dynes cm -2 , in a sample of galaxies with starburst-driven galactic superwinds. 
The derived P m in values for Mrk 6 are two to three orders of magnitude lower (see Table 3). 
If we assumed that the bubbles were powered by superwinds, then this low V min could imply 
that there is a large thermal contribution to the pressure, or the equipartition condition does 
not hold, or the bubbles are not powered by typical superwinds. In their work on the X- 
ray binary SS433, Blundell et al. (2001) have suggested that powerful accretion disk winds, 
which arise perpendicular to the radio jet, can also contribute to the radio luminosity. An 
accretion disk outflow cannot be ruled out in Mrk 6. 

In summary, the superwind model has the advantage that the complex geometry of the 
radio emission and the different features can be understood as being determined by their own 
local environments. The disadvantage of the model is that it cannot explain the similarities 
between the two pairs of the bubbles, if the similarities are not simple coincidences. However, 
the most important concern for this model stems from the derived energetics - the upper 
limit to the SFR provided by Spitzer observations indicates that a starburst alone cannot 
meet the energy requirements for the creation of the bubbles in Mrk 6. 
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4-2.2. AGN-driven Bubble Scenario 

The distorted appearance of many Seyfert jets has lead to the suggestion that the jet 
diverts due to interaction with molecular gas, within the central hundred parsecs (e.g., Gal- 
limore et al. 1996; Bicknell et al. 1998; Gallimore et al. 2006). The jet eventually terminates 
in a bow shock which plows into the surrounding ISM, resulting in bow-shock-like structures 
on kiloparsec scales, similar to the structures observed in Mrk 6. Assuming that the effi- 
ciency (e) with which the total jet energy is tapped to produce radio luminosity is 1% (see 
for e.g. O'Dea 1985), the jet with a radio luminosity of ~ 10 40 ergs sec -1 (see Table 3), and 
a kinetic luminosity (^^) of ~ 10 42 ergs sec -1 , would require ~ 3 x 10 5 years to deposit 
all its energy into the nuclear ISM and produce the derived minimum energy of ~ 10 55 ergs 
for the outer bubble. The jet would require only ~ 7 x 10 4 years to produce the derived 
minimum energy for the inner bubble. Although this model seems energetically feasible it 
does not naturally explain the two edge-brightened nested radio structures or the fact that 
the inner jet is perpendicular to the inner set of bubbles. 

We note that the pressure- volume energy of the outer bubble is around a hundred times 
larger than the derived minimum energy. Therefore, if the equipartition condition holds, 
the energy contribution from the thermal particles to the bubbles, must be a hundred times 
greater than from the relativistic electrons. 

AGN-driven plasmons have also been proposed to explain bubble-like structures in 
Seyfert galaxies (Pedlar et al. 1985). Models of radio emission from idealized plasmons 
(e.g., Shklovskii 1960) suffer severe adiabatic losses resulting in their radio luminosity de- 
creasing as r~ 5 , where r is the radius of the plasmon. Assuming adiabatic losses in a spher- 
ically expanding radio plasmon of relativistic gas, the final luminosity and internal energy 
of the plasmon is given by L = L (r /r) 5 and E = E (r /r), respectively, where L and E 
are the initial luminosity and internal energy of the plasmon with initial radius tq. Using 
L = 7.4 x 10 38 ergs sec -1 and E = 5 x 10 53 ergs at a radius r =7. 5/2=3. 75 kpc for the outer 
bubble (see Table 3 for <fi = 10 -3 ), we derive a large initial luminosity and internal energy of 
L ~ 6 x 10 56 ergs sec" 1 and E ~ 2 x 10 57 ergs, at the radius r =l pc. Similarly for the inner 
bubble with r=3.8/2=1.9 kpc, we derive L ~ 2 x 10 55 ergs sec" 1 and £ ~ 2 x 10 56 ergs, at 
r =l pc. However, typical bolometric luminosities for Seyferts galaxies fall in the range of 
10 45 — 10 46 ergs sec -1 (Padovani & Rafanelli 1988), making the L estimate unreasonable. 
Alternatively, an energy equivalent of 10 5 years of bolometric luminosity must be released at 
the radio frequency alone, at one instant. A single plasmon model for the bubbles of Mrk 6 
therefore seems to be energetically implausible. 
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4-2.3. Precessing Jet Scenario 

In this model, an episodic precessing jet is the prime mover that has produced all of 
the radio structure visible in Mrk 6. Assuming a jet velocity of ~ 10 4 km s" 1 (Bicknell 
et al. 1998; Ulvestad et al. 1999), and the bubbles-edges to be a precessing jet, the outer 
bubbles would have resulted in ^ 7 x 10 5 years and the inner bubbles would need ^ 1 x 10 5 
years to form. The gap between the size-scales of the two pair of bubbles would imply 
different episodes of ejection separated by ^ 6 x 10 5 years. Therefore, the outer bubbles 
would be the relic of the radio jet emission from about 10 6 years ago as it precessed in 
the northeast-southwest direction, probably following an accretion event. The precessing 
jet subsequently changed direction and pointed in the east-west direction about 10 5 years 
ago, possibly following another accretion/ejection episode. If the jet velocity is closer to 
10 3 km s _1 (e.g., Whittle et al. 2004), the time-scales for the formation of the outer and 
inner "bubbles" would be ^ 7 x 10 6 years and ) 1 x 10 6 years, respectively. 

With this picture in mind, we derived the proper motions of two precessing jets with 
j3 = 0.03, corresponding to the two phases of activity of the AGN, using Eqns. 1-4 of 
Hjellming & Johnston (1981). In Figure 6, we have plotted the proper motions of the 
precessing jet projected on the sky (i.e., the /i Q and {is parameters as described in Hjellming 
& Johnston, 1981), for two different jet inclinations and precession cone opening angles. 
The parameters pertaining to this model are listed in Table 4. The model suggests that 
the northern outer and the eastern inner jets are approaching jets (denoted by solid lines in 
Fig. 6), while the southern outer and the western inner jets are receding jets (denoted by 
dashed lines in Fig. 6). Further, the outer jets precess in the clockwise sense, while the inner 
jets precess in the counterclockwise sense. The jet flips its axis by ~ 140 degrees between 
the two epochs. 

The similarity between the precessing jet model and the radio images is striking. In 
the model, assuming only adiabatic losses and minimal radiative losses in the radio emitting 
plasma, the right edge of the northern outer "bubble" and the left edge of the southern outer 
"bubble" would both be brighter and have a flatter spectral index, as it would correspond to 
the most recent position of the jet. Elsewhere the spectra would have steepened due to the 
loss of energy in the synchrotron-emitting particles. The same would apply to the bottom 
edge of the eastern inner "bubble" and the top edge of the western outer "bubble" . As we see 
in Fig. 6 and Fig. 4, this is in good agreement with the observations. Further, the precession 
model easily reproduces the elongated morphology of the outer bubbles and the north-south 
S-shaped MERLIN jet - it is the projected terminal position of the jet that formed the inner 
"bubble" structure. 

The precessing jet model yields precession periods of 3 x 10 6 and 4 x 10 5 years for the 
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outer and inner jets, respectively (Table 4). Further, the jets were "on" for approximately 
one precession period in each episode of activity. The similarity in the precession and the 
"activity" period is perhaps suggestive of a link between the two phenomena - the accretion 
event itself could be the cause of the precession. It is very interesting that Natarajan & 
Pringle (1998) find the timescale for the accretion disk/blackhole alignment for a 10 6 — 10 7 M Q 
blackhole (see for example Peterson et al. 2004, for Seyfert blackhole estimates) accreting at 
one-tenth of the Eddington limit, to be of the order of 6 — 7 x 10 5 years. 
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Fig. 6. — A comparison of the morphology of an episodically-powered precessing jet (top right) 
and the 20 cm radio image (top left). (Right) Plots of the proper motions of the precessing jet 
projected on the sky (i.e., the (i a and parameters as described in Hjellming & Johnston, 1981). 
The axes are in arcseconds, the solid line represents the approaching jet, while the dashed line 
represents the receding jet. The outer and inner "bubbles" are jets from two different phases of 
activity of the AGN. Jl, CJ1, J2 and CJ2 indicate the jet and counterjet from the two phases 
and their approximate end-positions. In the precession model, the right edge of the northern outer 
bubble and the left edge of the southern outer bubble (both marked "New") would be brighter 
and have a flatter spectral index, as it is corresponds to the most recent position of the jet. The 
same applies to the bottom edge of the eastern inner bubble and the top edge of the western outer 
bubble (bottom panel). This agrees with the image on the left. Finally, the S-shaped MERLIN jet 
is easily reproduced in this model. 
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Table 4: Precessing jet parameters 



Structure 


i 




X 


P 


P 


t 




(cleg) 


(cleg) 


(cleg) 




(yr) 


(yr) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Outer jet 


153 


24 


-30 


0.03 


3xl0 6 


1.7xl0 6 


Inner jet 


9 


9 


-30 


0.03 


4xl0 5 


4xl0 5 



Note. — Col.l: In the precessing jet model, the outer and inner jets correspond to the the outer and inner 
"bubbles" in Mrk 6, Col. 2: Inclination angle of the jet axis, Col. 3: Half-opening angle of the precession cone, 
Col. 4: Angle needed to rotate the geometrical model so that the axes coincide with the true north and east 
(see Hjellming & Johnston 1981), Col. 5: Jet speed with respect to speed of light, Col. 6: Precession period 
in years, Col. 7: Total time the jet was "on" in years, assuming the jet speed given in Col. 5. The jet sign 
parameter, Sj e t — +1 was used for the approaching northern jet and Sj e t = —1 for the receding southern 
jet. The rotation sign parameter, s rot = —1 was used to imply clockwise motion of the precessing outer jet 
while s rot = +1 was used for the counterclockwise precessing inner jet. 
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In order to test the precession model, we derived a spectral age using the difference in 
the spectral index values taken from the two opposite edges of the eastern inner "bubble" 
(marked 'New' and 'Old' in Fig. 6, bottom right), under the assumption that the spectral 
steepening is a result of ageing. Using the Jaffe & Perola (1973) model for spectral ageing 
and following the analysis described in Myers & Spangler (1985), we derived a radiative 
age of 4 x 10 5 years, assuming an equipartition magnetic field of ~ 90/iG (Table 3). This 
estimate is a good match to the timescale the inner jet was "on" in the precessing jet model, 
providing further support to the model. 

The fractional polarization values along the edges of the bubbles range from 20% to 40%, 
values that are typically observed in AGN jets. Typical values for the measured polarization 
of supernova remnants, on the other hand, are less than 10% on all scales (eg., Gaensler & 
Wallace 2003; Lazendic et al. 2004; McConnell et al. 2006). The polarized emission indicates 
that the southern side of the outer and the western side of the inner bubbles is more polarized 
than the opposite sides. However, the radio polarization seems to be anticorrelated with the 
presence of the emission-line gas (see Kukula et al. 1996), suggestive of depolarization of the 
northern and eastern bubble by ionized gas. Such an anticorrelation has also been observed 
in some radio galaxies (Heckman et al. 1982; Liu & Pooley 1991). We have estimated V min 
for volume elements along the edges of the bubbles, assuming a volume filling factor of 
unity (see Table 3). In the precessing jet scenario, these correspond to sections of the radio 
synchrotron jet. We find the minimum pressure to be ~ 3 — 6 x 10~ n dynes cm -2 and the 
minimum pressure magnetic field to be ~ 16 — 26 /zG for sections of the inner bubble edge. 
The minimum pressure and magnetic field strengths are lower by an order of magnitude in 
portions of the outer bubble edge. We find that these values are similar to those obtained 
for radio jets in radio galaxies (Bridle et al. 1981; Bridle & Perley 1984) and Seyferts (e.g., 
Pedlar et al. 1985). 

The 6 cm MERLIN image of Mrk 6 reveals a gentle S-shaped knotty radio source (easily 
reproduced in our precession model, see Fig. 6) which has often been suggested to be the 
signature of a precessing jet. Episodic ejection of material into a two-sided, precessing 
(or wobbling), collimated outflow or jet has been proposed to explain the point-symmetric 
morphology of planetary nebulae (Livio & Pringle 1997; Raga et al. 1993; Lopez et al. 1993; 
Soker & Livio 1994; Livio & Pringle 1996). Specific examples include the pre-planetary 
nebulae IRAS 16342-3814 (Sahai et al. 2005) and Hen 3-1475 (Velazquez et al. 2004), the 
planetary nebulae K3-35 (Miranda et al. 2001) and M2-9 (Schwarz & Monteiro 2004), the 
AGB star W43A (Imai et al. 2002), the symbiotic star CH Cygni (Crocker et al. 2002) and 
Herbig-Haro type jets (Lim 2001). 

The S-shaped symmetry in jets has been observed in Seyfert galaxies like Mrk 3, 
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NGC 4151, NGC 5256, NGC 2685 and radio galaxies like 3C 305, 3C 345, 3C 294, 3C 334, 
and 3C 120 (Kukula et al. 1993; Pedlar et al. 1993; Heckman et al. 1982; Caproni & Abra- 
ham 2004b; Caproni et al. 2004). Most Seyfert galaxies observed on arcsecond-scales, do not 
show extended jets (e.g. Nagar et al. 1999). However when jet-like features are observed, 
they more often than not, exhibit a curved geometry (e.g. Ulvestad & Wilson 1984b; Kukula 
et al. 1995; Nagar et al. 1999; Thean et al. 2000). Indeed, 8 or S-shaped radio structures, 
sometimes along with large misalignments between radio structures on different spatial scales 
have been observed in the Seyfert galaxies NGC 2992, NGC 3079, NGC 5548, NGC 4051, 
NGC 6764, NGC 2110, NGC 4151, NGC 7469 (Wehrle & Morris 1988; Duric & Seaquist 
1988; Baum et al. 1993; Nagar et al. 1999; Thean et al. 2001). 

Further support for the precessing jet model comes from the emission line images ob- 
tained with the Hubble Space Telescope (Capetti et al. 1995), where the [OIII] and [Oil] 
emission of Mrk 6 follows the curve of the S-shaped jet far into the southern edge of the 
inner bubble. It would be difficult to explain this assymmetry in emission-line morphology 
with an expanding shell of starburst superwinds. The entire narrow-line region itself extends 
much beyond the extent of the radio source. A similar curvature in the narrow-line region 
and an ENLR much larger than the radio source have been observed in the Seyfert 2 galaxy, 
Mrk 3 (Kukula et al. 1993) and the radio galaxy 3C 305 (Heckman et al. 1982). Similar to 
the suggestions made for these galaxies (see Kukula et al. 1993), we interpret that the small- 
scale emission-line region is powered by the radio jets while the ENLR may be ionized by the 
radiation from the nucleus. However, in keeping with the picture of an episodic precessing 
jet, it is possible that some of the emission in the ENLR is a result of an earlier interaction 
with the radio jet. 

Therefore, based on all the above findings, we believe that the complex radio morphology 
of Mrk 6 can be understood in terms of an episodically-powered precessing jet. Mrk 6 may 
therefore be a unique object only in the sense that we are able to observe the source in its 
short-lived phase of episodic activity and before its radio emission from the previous epoch 
has faded away. Hence, similar radio structures may be found in other Seyfert galaxies going 
through such a phase, and when they are observed with sufficient resolution and sensitivity 
to image the relic emission from an earlier phase of activity. 



4-2.4- Implications of an Episodic Precessing Jet 

An episodically-powered precessing jet model in Seyfert galaxies implies that radio activ- 
ity is a short-lived phenomenon in the lifetime of a Seyfert galaxy. Further, multiple episodes 
of ejection are possible. These might be caused by independent accretion events. Based on 
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the fraction of spiral galaxies which are Seyferts, Sanders (1984) estimated the statistical 
Seyfert lifetime to be of the order of 3 — 7 x 10 8 yrs. Further, based on the physical extents 
of the narrow emission line regions and radio jets, Sanders concluded that Seyfert nuclear 
activity must be short-lived stochastic accretion events with no single episode lasting longer 
than 10 6 yrs. Although the estimate of the fraction of Seyferts in spirals has increased from 
a few percent to ~ 10% following more recent surveys (Ho et al. 1997), thereby increasing 
the statistical lifetime to 10 9 yrs, the episode lifetime remains the same. Random jet orien- 
tations in Seyferts occur, in his picture, because each accretion event is accompanied by its 
own independent accretion disk in no preferred plane. Based on the NICMOS imaging of 
a large sample of Seyfert galaxies, Hunt & Malkan (2004) also reached the conclusion that 
Seyfert activity may be prompted by a disturbance. 

The precessing jet model further implies that the jet ejection axis is easily perturbed in 
Seyferts. An accretion disk which is irradiated by a central AGN can become unstable to 
becoming warped, resulting in the wobbling of the jet (Pringle 1996, 1997; Livio & Pringle 
1997) and severe misalignment of the jet axis and the axis of the outer accretion disk. Apart 
from radiative feedback that may warp a rotating accretion disk, jet precession could arise 
due to near-Eddington accretion events that may alter the spin of the supermassive black 
hole, in turn affecting the orientation of the radio jet (Rees 1978; Scheuer & Feiler 1996; 
Natarajan & Pringle 1998); binary blackholes (Caproni & Abraham 2004a) and blackhole 
mergers, that might result in short timescale redirection of the jet axis (Merritt & Ekers 
2002; Merritt & Milosavljevic 2005). Possibly due to the smaller blackhole masses in Seyfert 
galaxies compared to typical radio galaxies (e.g., McLure & Dunlop 2001), the spinning 
blackhole axis and consequently the radio jet axis is perturbed easily when accretion events 
take place, resulting in precessing jets and switches in jet-orientation. 

It must be noted that the host galaxy of Mrk 6 shows no conspicuous signatures of 
galaxy-galaxy interaction or merger-activity, which could lead to a binary black hole in the 
centre (e.g., Begelman et al. 1980). In fact, Mrk 6 seems to be located in a sparse environment 
(Dahari 1985; Pfefferkorn et al. 2001), although the possibility of a major merger more than a 
dynamical timescale (~ 10 9 yrs) ago, cannot be ruled out. Note that the emission-line region 
in Mrk 6 is extended along more than one axis (Fig. 2). However, as Dennett-Thorpe et al. 
(2002) point out from their study on 'winged' radio galaxies, there remain two possibilities 
that could result in jet realignment without a second black hole. (1) An ingested dwarf galaxy 
may not leave observable traces that a larger merger would, and yet provide adequate mass 
and angular momentum to cause jet re-orientation. (2) The jet axis may not be determined 
by the black hole alignment, but be strongly influenced by the accretion disk. Instabilities 
in the inner disk could then cause switching of the jet direction over short timescales. 
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A consequence of the changes in jet orientation is that the jet may interact with the 
molecular torus, stirring up the constituent matter. This would be consistent with the 
variable absorbing gas column density inferred through XMM-Newton X-ray spectroscopic 
observations of the central regions of Mrk 6 by Schurch et al. (2006). 

5. SUMMARY AND CONCLUSIONS 

We have carried out an extensive radio study with the VLA on the Seyfert 1.5 galaxy 
Mrk 6 and imaged a spectacular radio structure in the source. The radio emission occurs on 
three different spatial scales, from ~7.5 kpc bubbles to ~1.5 kpc bubbles and a ~1 kpc radio 
jet, all lying roughly orthogonal to each other. All the radio structures appear to originate 
from near the AGN core. 

We have considered three models to explain the complex radio morphology of Mrk 6. 
In the first model, we have suggested that the bubble-like radio structures are the result of 
superwinds ejected by a nuclear starburst. However, recent Spitzer observations of Mrk 6 
provide an upper limit to the star formation rate of ~ 5.5 M© yr~ x , an estimate lower by 
an order of magnitude than the SFR of ~ 33 M Q yr -1 derived assuming that the bubbles 
are a result of starburst winds energized by supernovae explosions. Thus, a starburst alone 
cannot meet the energy requirements for the creation of the bubbles in Mrk 6. 

We have also considered a model wherein the bubbles are a result of diverted AGN- 
driven jets. The small-scale radio jet gets diverted and decelerated through its interaction 
with the surrounding interstellar medium, thereby depositing all its energy into the ISM 
in about 10 5 years. The energetics however, do not favour the AGN-driven single plasmon 
model proposed by Pedlar et al. (1985). 

In the third model, we attribute the complex radio structure to be the result of an 
episodically-powered precessing jet that changes orientation. The attraction of this model 
is that it can naturally explain the complex radio morphology. It is also consistent with the 
energetics, the spectral index, and the polarization structure. Recent XMM-Newton X-ray 
observations of a variable absorbing gas column density could be inferred as in-keeping- 
with the changing jet orientation in Mrk 6 which causes the jet to interact with the torus, 
thereby stirring up the gas/dust. Radio emission in this scenario is short-lived phenomenon 
in the lifetime of a Seyfert galaxy which results due to an accretion event. Further, possibly 
due to the smaller blackhole masses in Seyferts compared to typical radio galaxies, the 
spinning blackhole axis and consequently the radio jet axis is perturbed easily when accretion 
events take place, resulting in precessing jets and switches in orientation. Multi-wavelength 
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comprehensive radio studies of Seyfert galaxies would therefore reveal many more sources 
with similar complex radio structures as Mrk 6. 
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